Quasi Periodic Oscillations in X-ray Neutron Star 
and Probes of Perihelion Precession of General Relativity 
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We ascribe the twin kilohertz Quasi Periodic Oscillations 
(kHz QPOs) of X-ray spectra of Low Mass X-Ray Binaries 
(LMXBs) to the pseudo-Newtonian Keplerian frequency and 
the apogee and perigee precession frequency of the same mat- 
ter in the inner disk, and ascribe 15 - 60 Hz QPO (HBO) to 
the apogee (or perigee) precession and its second harmonic 
frequency to both apogee and perigee precession in the outer 
disk boundary of the neutron star (NS) magnetosphere. The 
radii of the inner and outer disks are correlated each other by 
a factor of two is assumed. The obtained conclusions include: 
all QPO frequencies increase and frequency difference of twin 
kHz QPOs decreases with increasing the accretion rate. The 
obtained theoretical relations between HBO frequency and 
twin kHz QPOs are simlilar to the measured empirical for- 
mula. Further, the theoretical formula to calculate the NS 
mass by the twin kHz QPOs is proposed, and the resultant 
values are in the range of 1.4 to 1.8 Mq. QPOs from LMXBs 
hkely provide an accurate laboratory for a strong gravitational 
field, by which a new method to determine the NS masses of 
LMXBs is suggested. 

PACS numbers: 97.60.Jd, 97.80.Jp, 04.80.Cc 



Until now about twenty LMXBs have been discover- 
ied to exhibit kHz QPOs, include one X-ray pulsar, with 
the Rossi X-ray Timing Explorer (RXTE) since the early 
1996, briefly after its launch The Z sources (Atoll 

sources), which are high (less) luminous neutron-star low- 
mass X-ray binaries [Q , typically show four distinct types 
of QPOs [Q. At present, these are the normal branch 
oscillation (NBO) ~ 5 — 20 Hz, the horizontal branch os- 
cillation (HBO) z/ff ~ 15 - 60 Hz §, and the kHz QPOs 
iy2{t^i) — 200 — 1200 Hz that typically occur in pairs, 
upper frequency 1^2 and lower frequency vi. In several 
Atoll sources, nearly coherent i^turst — 330 — 590 Hz os- 
cillations have also been detected during thermonuclear 
Type I X-ray bursts, which are considered as the spin 
frequency of NS or twice of them. All of these QPOs 
but the burst oscillation have centroid frequencies that 
increase with inferred mass accretion rate M, and are 
tightly correlated with each other (5|J|]. In some cases, 
these correlations appear to depend weakly on the other 
properties of the sources |^,^. For example, the frequen- 



cies 1^2 and h'l, as well as the frequencies 1^2 and fol- 
low very similar relations in five Z sources |^,^. How- 
ever the frequency separation between the upper and the 
lower kHz QPOs Av = V2 — vi decreases systematically 
with instantaneous M in some cases, e.g. Sco X-1 [0,0, 
4U1608-52 and 4U1728-34 [|l],||], then in the latter the 
observed coherent burst frequency 364 Hz is higher than 
its maximum Av ~ 355 Hz |^,^ . 

A number of theoretical models have been proposed 
to account for the QPO phenomena in X-ray NS sys- 
tems. For the high frequency of kHz QPOs and its pro- 
portional relation to the accretion rate, simply, the upper 
kHz QPO V2 is considered to originate from the Keple- 
rian orbital frequencies at the preferred radius close to 
the compact object, which exhibit the inner accretion 
flows, however the lower kHz QPO vi is attributed to 
the beat of such frequency with the stellar spin Vs . 
Recently, the general relativistic effects are paid much 
attetion to account for kHz QPOs pT| , p^ , which can ex- 
plain the varied kHz QPOs separation Av. Then the 
varied Au seems to exclude the simple beat explanation 
for kHz QPOs [Qj^], which predicts a constant frequency 
separation. However, it is also a common phenomenon 
for the separation Av to be approximately constant. The 
range of Av is also quite narrow across different sources 
(Ai/ w 250 — 360 Hz) with a nearly coherent frequency 
of Vburst ~ 330 — 590 Hz. In some cases (such as 4U 
1702-43 and KS 1731-260) the burst frequency is consis- 
tent, to within the errors, with the frequency separation 
Ai/, or twice its value 2Ai/. So the approximately con- 
stant separation in some sources seems also to make the 
purely GR precession model in difficulty. Although many 
other viable new ideas are also proposed (l^Jl^ , there has 
not yet been one model satisfactorily to explain all ob- 
servational QPO phenomena of LMXBs until now. The 
mechanisms of kHz QPOs of LMXBs are still debated 
and open problems. 

Nonetheless, HBO frequency [vh ~ 15 — 60Hz), first 
discovered in GX 5-1 in 1985, is interpreted to be the 
beat frequency between the Keperian frequency of the 
disk and the stellar spin frequency by the standard beat 
frequency model (BFM) ^ . 

It is commonly accepted that the kHz QPOs produce 



close to the innermost stable orbit (IMS) or the surface of 
the NS, which will provide a probe to detect the accretion 
flow around the non-Newtonian strong gravity region. 

Here we concentrate on the explanation of kHz QPOs 
and HBOs (for Atoll sources, 15 - 60 Hz QPO is sup- 
posed to be the same mechanism as HBO of Z sources 
1^ ) , and neglect the detail of the physical mechanism for 
QPO production. We ascribe 1^2 to the so-called pseudo- 
Newtonian Keperian frequency proposed by Paczyhsky 
and Wiita |l|,0 V2 ^ i^K = {GM/ATT'^r^Y''^{l - 
2GM/rc^)~^ , which is usually applied in the black hole 
accretion disk and coalecing of compact objects p^ ] 
to mimic the general relativistic effect, and to the 
apogee -|- perigee precession frequency ( ^^^2*^ rad per 
revolution) depicted by the post-Newtonian scheme |1^] 
1^1 — ^^¥-1^2- We suppose that both V2 and origi- 
nate from the same matter in the inner disk boundary. 
Further, contrary to the standard description to account 
HBO for BFM ||l^, as well as for the nodal precession 
[0 of Lense-Thirring effect in the inner disk, we assume 
that VH is an apogee (or a perigee) precession frequency 
in the outer disk boundary (magnetosphere boundary). 
Then, on the reasons for deliminatimg BFM to account 
for HBO, we have the following arguements. First, the 
simple BFM cannot predict the second harmonic fre- 
quency iyH2, which was detected in Sco X-1 with 90 Hz 
corresponding to HBO 45 Hz Q. Second, vh, V2 and 
vi correlations seems to weakly depend on the stellar 
spin frequency 0,^. Third, BFM should predict the null 
or extremely low vh if the accretion rate is extremely 
low, where the inverse of the spin-up torque should arise 
vh to be inversely correlated to the accretion rate. But 
this phenomenon has never been detected in both Z and 
Atoll sources, and the lowest detected HBO vh ~ 15Hz 
with the lower accretion rate hinted in color-color dia- 
gram of X-ray spectra Well, here the proposed pe- 
riastron precession model (PPM) can ascribe the second 
harmonic i'h2 to the precession frequency of both apogee 
and perigee at the outer disk boundary, which is twice of 
HBO {iyH2 = '2vh). 

If this PPM interpretation were confirmed, QPOs 
would provide a probe into the test for strong field GR, 
by which we can determine or constrain the NS parame- 
ters, such as mass, radius and magnetosphere radius, as 
well as the magnetic field strength. We set c = G = 1 
throughout this paper. 

On account of the complexity of the accretion flow 
close to the innermost stable orbit Ri — 3i?s, where 
Rs — 2GM/c^ — 3m (km) is the Schwarzschild radius 
calculated with the gravitational mass M (m is the mass 
in unit of solar mass), the motion of the disk matter will 
be influenced by not only gravitaional field but the star 
magnetic field, which should be complicated. For sim- 
plicity, we assume that the disk matter to exhibit QPO 
are mainly dominated by the Schwarzschild gravitational 
field with a slightly eccentric orbit (e ~ 0). Therefore, 



these QPO frequencies are conveniently arranged as fol- 
lows with the consideration of the possible modification. 
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where r is the inner disk radius, and Rco is the outer disk 
radius where it is corotating with the magnetosphere. 
i/fem is the pseudo-Newtonian Keplerian frequency of 
the outer disk, ai, 012, a^, ctfem 1 are the modifica- 
tion coefficients on account of the orbit eccentricity and 
other approximated uncertainty, which might be almost 
unity under the ideal condition with the circular orbit. 
^i(fe) = 0,1,2 represent the possible beat modes be- 
tween kHz QPOs and HBO. (/) is a scaling parameter to 
connect the inner and outer disks, but here we suppose it 
to be 0.5 for the reason of the best fitting, y is the ratio of 
the Schwarzshild radius to the instantaneous inner disk 
radius, wich is proportionally related to the fractal accre- 
tion rate, Mr = -r^, the ratio between the instantaneous 

' (M) ' 

accretion rate M and the average accretion rate (M), 
and the latter determines the corotation magnetosphere 
radius i?co = 1.9 x 106(cTO)Bg/^(M)-2/7mi/7^i2/7 j-^^^j^ 

the nearly unity fastness [ p^|p^ ) , where and Rq are 
star magnetic field and radius in unit of 10* G and 10^ 
cm respectively. 

The relations vh vs. V2 and l^u vs. ui are plotted in 
Fig.l and Fig. 2, together with the well measured source 
samples, and it is shown that the agreement bewteen the 
model and the observed QPO data is quite well for the se- 
lected values of the NS mass from 1.4 to 1.8 M©, which is 
the solely free parameter in the Eqs.(l-5). From Eqs.(|]), 
(^) and (^, we can derive the thoeretical relations be- 
tween QPO frequencies in the following if we set the ideal 
parameter conditions with ai = a2 = ah ~ ctkm ~ 1 and 
kx=k2 = 0. 



^ 44.2 (i/z) (^)[l-0.1(f|j-)2/5] , 



(6) 



51.5 {Hz) mV3(^)5/3[l _ i(^)2/3]2/3 ^ (7) 



^1 ^ 595 {Hz) m2/3(^)5/3[l _ l{^fl^fl^ . (g) 

Further, we can also obtain the NS mass formula repre- 
sented by the QPO frequencies without consideration of 
the other modification parameters 



m = 2.17x(^)3/2(^)-i[l-i^]-i, (9) 

m 1.56 X (^)3/2(_^)-5/2[i _ 6^]-l ^ (10) 

For Z sources (almost Eddington accretion) , — 1 cor- 
responds to the maximum QPO frequencies, then we can 
determine ymax by the observed iy2max and calculate the 
corotation magnetosphere radius, as well as the spin fre- 
quency. It is stressed that the observed maximum upper 
kHz QPOs are most likely the maximum kHz QPO fre- 
quencies for Z sources but they are not neccessary. How- 
ever for Atoll sources (lower accretion rate) , they have no 
above properties shared by Z sources. But, we can apply 
the measured h'2max to predict the possibly upper limit 
for the burst frequency (spin) of Atoll sources. Moreover, 
for the known burst frequency (spin) , the corotation Ke- 
plerian frequency corresponds to = 1, by which we 
can determine the magnetosphere radius and magnetic 
field strength of the source. 

The estimation of the NS magnetic field strength 
III is given hy B = {^f'^Bf with Bf = 4.3 x 

10« (G) {{M)IMEddf'^m^/^R~^'^, where R is star ra- 
dius. The calculated NS parameters of LMXBs are list in 
TABLE I for six Z sources and eleven Atoll sources. The 
obtained magnetic fields for both Z and Atoll sources are 
about 2 X 10^ G and 2 x 10* G respectively, which are 
consistent with the originally suggested values |^,^. It 
is also interested that the derived spin frequencies are in 
the range of 289 (Hz) - 356 (Hz) with hinted magneto- 
sphere radius in the range of 37 (km) - 47 (km), and the 
homogenous span of spins of LMXBs is consistent with 
the period distribution of millisecond pulsars. However 
the reason for this homogeneity is still in studying [ p2| . 

Further, we can define the maximum kHz QPO fre- 
quency (fci — k2 — Q mode) at the innermost stable orbit 

Rl = 3Rs or y=l/3, then V2max = Vlmax = V2{v = 

1/3) — 1625 X (2/m)(iJz), which is much higher than 
the measured maximum kHz QPO frequency 1228 Hz 
(4U1636-53). The higher V2max seems to hint the satura- 
tion of V2 will not happen below 1200 Hz if the levelling 
off of the X-ray signal orginates from the IMS. The re- 
cent observation on the kHz QPO saturation obtained 
the debated conclusions [pH23|. 

It is easy to see in FIG. 2 that the kHz QPO separa- 
tion Az^ vs. vi for two well measured sources Sco X-1 
and 4U1728-34 ||l|,0jD is plotted, together with the pro- 
posed PPM theoretical curves. It is apparent that for NS 
masses at about 1.45M0 and 1.7OM0, the simple model 
outlined above is in qualitative agreement with the mea- 
sured values, including the decrease of for increasing 
vi. For the separation /S.v to be approximately constant, 
we argue that the beat mode with fci = 1 and fc2 = 
can mimic the almost constant separation in the range 
of vi = 500 ~ 900(i/z), corresponding to the parameter 
range oiv = 0.2 ~ 0.25. For the mode {k-\ =1, ki = 0), 



the separation becomes ^v{y) ~ (1 — 2.75i/)j^i(fci = 0), 
and we will find Az^(y = 0.2)/Aj/(2/ = 0.25) ~ 1. If the 
spin frequency is defined by the corotation frequency 
of magnetosphere-disk of about 40 (km) with y=0.25, 
i.e. Us — Vkm{y — 0.25), then we obtain that the ra- 
tios Ai/(y = Q.T)/vs ^ 1.01 and Au{y = 0.25)/u, ~ 0.99, 
both ratios are independent of the NS mass m ! This may 
be the reason why for a few sources the measured Au are 
similar to the burst frequencies within the errors. This 
coincidence refiects the homogenous magnetospheres of 
LMXBs in the range of about 40 (km). 



TABLE I. The parameters of six Z sources 
and eleven Atoll sources 



Source 


'^27nax m Rco Vs B^ 

(Hz) (Mo) (km) (Hz) (G) 


GX 5-1 


930^^ 1.71 46.7 289 22 


GX 340+0 


945° 1.43 42.5 297 18 


Cyg X-2 


1005 1.50 42.1 313 17 


GX 349+2 


1020 1.78 45.6 314 21 


GX 17+2 


1085 1.63 42.3 334 18 


Sco X-1 


1130 1.75 43.0 345 19 


4U0614+09 


1145 1.41 38.2 356^ 1.5 


4U1608-52 


1090 1.65 42.4 335^ 1.8 


4U1705-44 


1075 1.61 42.2 332^ 1.8 


4U1735-44 


1150 1.62 41.0 352^ 1.7 


4U1820-30 


1065 1.55 41.6 330^ 1.7 


4U1915-05 


1005 1.37 40.3 316^ 1.6 


XTEJ2123-058 


1140 1.70 42.2 349^ 1.8 


4U1636-53 


1.75 45.5 290* 2.1 


4U1702-43 


1.49 39.4 330* 1.6 


4U1728-34 


1.44 36.9 364* 1.4 


4U1731-260 


1.50 45.5 524/2* 2.0 



a : Inferred from the frequency separation of kHz 
QPOs. *: data from van der Klis Q, and references 
therein. ] : the possibly upper limit for the burst fre- 
quencies(or half of them). +: star radii are assumed 
to be 15 km with the Eddington accretion rate MEdd for 
Z sources and Q.QlMEdd for Atoll sources. 

In summary, we state that the model described here is 
simply and roughly one, and many physical details are 
neglected, such as NS spin induced gravitomagnetic ef- 
fect, NS quadrupole induced nodal precession ||ll[, the 
self-gravity of the disk, magnetosphere structure and 
magnetic axis inclination, the spiral-in effect of the ac- 
creted matter, the origination and inffuence of the non- 
zero eccentricity, etc.. Especially, the non-zero eccentric- 
ity should have somewhat effects on the QPOs, but the 
origination mechanism for this is still unknown. At least, 
we can speculate that the motion of the accreted matter 
in the disks might be not exactly described by a free test 
particle in a circular orbit of a purely gravitational field. 
Consideration of these factors will construct our future 
exploration and understanding of QPOs in LMXBs. 



Finally we stress that the proposed PPM model has 
also implications that the QPO phenomenon in LMXBs 
is likely to reflect the fundamental relativistic motion of 
matter in the vicinity of the NS strong gravitational field. 
As discussed in this paper, if the PPM were successfully 
to explain the observed QPOs in LMXBs, we not only for 
the second time uncover the GR periastron precession 
evidence in the NS besides that in the Mercury orbit 
of the solar system, predicted by Einstein, but also, by 
means of QPO, we developed a new method to determine 
the NS mass from Eq.(0) and Eq.(lO), as a contrast to 



determining NS mass in NS binary system such as PSR 
1913+16. 

The author is grateful to E. Ford, M. Mendez, R. 
Wijnands and W. Zhang for critical reading of the 
manuscript and helpful discussions. 
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FIG. 1. HBO frequency vh versus the upper kHz QPO 
frequency 1^2 for five Z sources of LMXBs. (cf. [4,5] and ref- 
erences therein). Error bars are not plotted for the sake of 
clarity. The model presents a well fitting for the nearly circu- 
lar orbit of NS mass from about 1.4 to 1.8 solar mass. 
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FIG. 2. kHz QPO frequency separation Ai/ versus the 
lower kHz QPO frequency U2 for Sco X-1 and 4U1728-34 
(cf. [1,6-8]). Error bars are not plotted for the sake of clar- 
ity. 



